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ABSTRACT 

An efficient diagonal scheme implemented in ail 
overset mesh framework has permitted the analysis of 
geometrically complex cavity flows via the Reynolds- 
averaged Navier-Stokes equations. Use of rapid hyper- 
bolic and algebraic grid methods has allowed simple 
specification of critical turbulent regions with an al- 
gebraic turbulence model. Comparisons between nu- 
merical and experimental results are made in two- 
dimensions for the problems of: a backward- facing step, 
a resonating cavity, and two quieted cavity configura- 
tions. In three-dimensions, the flow about three early 
concepts of the Stratospheric Observatory For Infrared 
Astronomy (SOFIA) are compared to wind-tunnel data. 
Shedding frequencies of resolved shear layer structures 
are compared against experiment, for the quieted cav- 
ities. The results demonstrate the progress of compu- 
tational assessment of configuration safety and perfor- 
mance 1 . 

NOMENCLATURE 


c 

speed of sound 

c p 

coefficient of pressure, 1 

Voo 

f 

frequency 

I\ 

ratio of convection by freest ream speed 

t 

mixing length 

L 

characteristic length 

w 

stage number 

m 

mass flow rate 

M 

Mach number 

P 

instantaneous static pressure 

PSD 

power spectral density, dB 


q velocity magnitude or dynamic pressure 

Q vector of dependent variables 


r 

velocity ratio, ^ L 

• 7 11 2 

Pc 

Reynolds number 

St 

SPL 

Strouhal number, — ^ — 

tii -r !i 2 

sound power level, (ID 

t 

time 

u , r, w 

Cartesian velocity components 

y, ~ 

Cartesian physical space coordinates 

a 

angle of attack 

9 

momentum t hickness 

X 

wavelength 

Pt 

eddy viscosity 

P 

density 

( 7 

spreading rate parameter 

n 

mean quantity 

< > 

root mean square quantity 

( r 

fluctuating quantity, / = / + f 

Subscripts 

a 

acoustic 

T 

V 

total quantity 
vor tical 

oo 

freestream quantity 


INTRODUCTION 

The effort reported here describes the progress of 
a computational approach for use in the design of a 
new airborne astronomical observatory. The existing 
Kuiper Airborne Observatory (KAO) has provided two 
decades of unique research capabilities. 1 Figure 1 de- 
picts the proposed successor, SOFIA, which will offer 
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ten times more resolution than the KAO. Assessment 
of the safety ami performance of this large cavity in a 
Boeing 7I7SP is the goal of this continuing elfort. To- 
wards this objective of providing design information, 
extensive (‘valuation of tin' numerical methods by com- 
parison against experiment is necessary. 

This report describes the status of the validation of 
the computational methods to date. In order to pro- 
vide a somewhat complete overview, both previously 
reported 2 and new results are incorporated herein. The 
two-dimensional cases discussed here include free shear 
layers, 3 a backward-facing step, 4 a resonating cav- 
ity, 5 and two quieted cavities. 14 The computed three- 
dimensional cases are compared to the wind tunnel data 
of Rose and Cooley.' The following sections address the 
method used to predict the unsteady Hows, including 
the grid generation and modelling of turbulence. 



Fig. 1: An early conceptualization of SOFIA 


METHOD 

The fluid field was computed via the Navier-Stokes 
equations using the diagonal scheme of Pulliam and 
Chaussee 8 implemented in the overset grid framework 
of Bcnek, Bulling, and Steger. 9 The equations were 
integrated through Euler implicit time marching and 
second-order spatial differencing with viscous wall con- 
ditions specified as no-slip, zero normal pressure gra- 
dient, and adiabatic. Information transfer across over- 
set mesh boundaries was implemented using trilinear 
interpolation of the dependent variable vector, Q = 
[pipu, pi\ pw,e] 7 . Computations were performed on 
the NAS Cray 2 and Y-MP, the CCF Cray Y-MP, and a 
Silicon Graphics 4D35TG workstation. The flow solver 
cost is 1 3// ,s /coll /stop on a Cray Y-MP. 

Turbulence Model 


win and I.omax. 10 ;is implemented by Bulling, 11 is de- 
scribed below using a flow in the (.r.i/) plane. 

The wall-bounded flows use the original Baldwin- 
Lomax model, with the addition of a variable F, uaT 
cutoff. The grids are chosen swell that a unique wall 
distance is readily available. 

The cavity aperture spanning shear layer uses an 
eddy viscosity developed using F(y) = y |-c|, as sug- 
gested by Baldwin and Lomax 10 for wake regions. This 
results in 


F tpd k e — C u'k 
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where specification of C w k is discussed below and the 
velocity difference is modified to be half the total veloc- 
ity difference between the streams in the specified shear 
layer region 

<1,1, f = v/(« 2 + - vA'' 2 + 


The free shear layer model is now given by 
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after dropping the Klebanoff intermit t-ency function. 
Unmodified model constants K = 0.0 1G8 and C cp = l.G 
are used. 

The magnitude of the eddy viscosity in the free shear 
layer model can be altered with C\ r k • Estimation of 
the proper value of C wk begins by using GortlcEs shear 
layer solution: 
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where u i and u 2 are the velocities of the slow and fast 
streams and £ is the similarity coordinate, as shown in 
Fig. 2. The spreading parameter <r is inversely related 
to the spreading rate, dbjdx , where 6 is a measure of the 
shear layer width. The value of the spreading parameter 
when the velocity of one of the streams is zero is cr {) . 

Gbrtlcr's solution can be used to determine the max- 
imum vorticity magnitude as follows: 
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In these computations, the slowly time- varying com- 
ponent of the flow is resolved, while rapid fluctuations 
are modelled. The algebraic turbulence model of Bald- 
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Now, using Prandtl's mixing length assumption and 
scaling laws for jet boundaries, eddy viscosity can also 
be expressed as 


fl t X pC 2 \u>\ 

,2 A " 

X ^ — 

= A'd/j/Ah 


( 2 ) 


where AN = tt-t- 
ibal 

Setting Eqs. 1 and 2 equal results in 


Cwk — (<7o A C cp v/tt) 


and only <7o remains to be specified. Empirical esti- 
mates of (7 0 range from 9.0 to 13. 5. 3, 12 13 For this series 
of cavity flow efforts <7o was set to 11.0, resulting in a 
value of C w k = 1-91. 


Grid Generation 

Computation of the loads generated by cavity flows 
requires accurate representation of the geometry as well 
as the flowfield. Typically, a significant effort, in grid 
generation is required before flow analysis can begin. 
Since matching zone faces arc not required for the over- 
set method used here, recent advances in algebraic 15 
and hyperbolic 14 methods can be used. Hyperbolic grid 
generation, which gives good spacing and orthogonality 
control, was used for the wall-bounded regions, while al- 
gebraic grids were used in shear flow regions including 
plumes and wakes. Advantages of this type of grid sys- 
tem include straightforward specification of the turbu- 
lent regions and allowance for independent refinement 
of each zone. This topology also permits the re-use of 
meshes for configuration studies. 


FREE SHEAR LAYER 

Numerical experiments were performed using a two- 
dimensional shear layer to determine sensitivities of 
mean and time-varying quantities to changes in time 
step size, fourth-order dissipation levels, and grid refine- 
ment. In addition, the algebraic turbulent shear layer 
model was partially validated through comparison with 
similarity solutions and experimental data. 

The computational domain for this case is shown to 
scale in Fig. 2a, where inviscid channel side walls and 
characteristic inflow and outflow boundary conditions 
were used. The boundary layers on the splitter plate 
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Fig. 2: Shear layer case: (a) velocity profiles of differing 
grid resolution compared to similarity and (b) variation of 
spread rate with velocity parameter 


and the shear layer were fully turbulent with a Reynolds 
number based on the mean velocity of the streams and 
the length of the splitter plate given as G.7 x 10 s . 

The results for three grid refinement levels are also 
shown in Fig. 2a along with Gortler's similarity solu- 
tion, where the velocity profiles were taken more than 
a thousand momentum thicknesses downstream of the 
trailing edge of the plate. The solution can be seen to 
become grid dependent when fewer than than 20 points 
span the layer. Eddy viscosity was observed to grow lin- 
early in accordance with the Clauscr formulation, while 
the solution w;is insensitive to fourth-order dissipation 
levels within the range 0.01 to 0.05. The Mach ratio for 
this case was 0.2/0. 8 with a = 20.7. 

Figure 2b compares the variation of spread rate with 
velocity ratio, and demonstrates that computed spread- 
ing rates are within the bounds of the experimental 
data. 

BACKWARD-FACING STEP 

A qualitative estimate of the errors to be expected 
in the three-dimensional flows can be gleaned from the 
backward-facing step problem. The computation was 
compared to the experimental data of Driver and Sceg- 
miller, 1 with test conditions matched at 

M ref = 0.128, He // = 3.44 x 10 l 

p r€ f = 1.10 k<j/m\ Prcf = 9.11 x 10 4 A 7 in 2 

Characteristic inflow and outflow boundaries were 
used, holding mass flow, total enthalpy, and flow angle 
fixed at the inlet while specifying pressure at the outlet. 
Both m and p rx u were iteratively adjusted to match two 
experimental conditions. First the overall pressure rise 
from x/H = -4 to x/H — 30 was matched to experi- 
ment, then m was modified to match the experimental 
momentum thickness at x/H = —4. Finally, viscous 
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adiabatic walls won' specified for doth the top and bot- 
tom of the channel. The unusual topology shown in 
Fig. 3 was used to replicate those required by the more 
complex configurations. 



Fig. 3: Backward-facing step: coarsened near field grids 

After dissipation of initial transients, the simulation, 
which was performed using a step size of 15//$, was av- 
eraged over 200 steps by increments of 20 steps. Using a 
particle convccting with the shear layer from separation 
to the reattachment point, then 200 steps corresponds 
to three characteristic time increments. The time mean 
field was sampled in 7//10 increments via bilinear in- 
terpolation of a valid donor cell to obtain the profiles 
shown in Fig. 4. 




Fig. 4: Backward-facing step: velocity profiles 

Agreement in maximum vorticity is seen in Fig. 4, 
however some discrepancies in profile shape exist. Us- 
ing time mean skin-friction levels, the predicted reat- 


tachment occurred at j*/// = G.3 while experimental 
data showed roat t aelmicnt at x/H = G.0. 

TWO-DIMENSIONAL RESONATING CAVITY 

In order to establish confidence in the numerical 
method, a two-dimensional cavity computation was un- 
dertaken to demonstrate and validate'* self-induced cav- 
ity resonance. 

The test conditions, specified to match experiment, 
were 

Mrn = 0.9, Rci = 1.47 x l() r> , L - 8 in. 

Poo = 0.40 kg/m 3 , p x = 2.9 x 10* X/m 2 

and a cavit} - length by depth, Lj D, of 2. Characteristic 
inflow and outflow conditions were specified along with 
a step size of Af — 1.97//$. Power spectra were com- 
puted using 8192 steps following the dissipation of the 
initial transients. 

Inspection of the computed cavity pressure history, 
shown in Fig. 5a, confirms the idealized feedback mech- 



Fig. 5: 2-d cavity: (a) pressure history and (b) power 
spectra comparison 

anism identified from Rossiter’s experiments. Briefly, 
the cycle begins with the propagation of an acoustic 
wave from the aft wall of the cavity to the forward bulk- 
head. Wave reflection from the forward wall causes the 
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shear layer to bow outwards. shedding vort icily. The 
d e {] ee t ed shear layer converts downs! ream and induces 
another cycle. This coupling of the acoustic and vor- 
tical fields is quantified hy Kossiter's empirical model, 
given in Fig. oh, which gives only feedback frequencies. 

In the frequency domain, comparison of Rossi t or *s 
data to present results indicate agreement in frequency 
at the peak magnitudes, as shown in Fig. ob. Mag- 
nitudes are higher for the present case by about 2 dB, 
which can be explained from dimensionality arguments. 
The solution was also found to be insensitive to second- 
order dissipation levels within the range 0.3 to 0.5. 
Figure 5b also shows the resonant modes predicted by 
Rossiter’s equation, showing that. A = 0.5G gives bet- 
ter prediction of the higher modes. Finally, the vertical 
knife edge schlieren images of Fig. 6 show the qualita- 
tive agreement between computed and observed 16 radi- 
ation patterns. 


Experiment, Karamched Present 

Fig. 6: Comparison of experimental and numerical 

schlieren images with knife edge vertical 

TWO-DIMENSIONAL TREATED CAVITY 

The effect of cavity geometry, particularly modifica- 
tion of the shear layer attachment region, is known to 
possess potential quieting capabilities. 17 The Army Air- 
borne Optical Adjunct (AOA), shown in Fig. 7, flight 
tested several passive and active quieting methods. 6, 18 
The purpose of the present numerical simulations is to 
determine if optical quieting methods, particularly aft 
ramp treatment and lip- blowing, could be accurately 
simulated. 

The grid cell size was specified as 0.83 in. in the 
st.reaimvise direction, chosen so that frequencies up to 
approximately 400 Hz would be resolved without signifi- 
cant numerical dissipation effects. 19 A time step of 44//s 
was fixed so that CFL « 1 in the streamwise direction 
within the shear layer. The numerical test, conditions 
were mat died to flight data: 

A/oo = 0.77, Rc l = 5.00 x I0 6 

= 0.2G2 kg/m 2 , Poo = 1.G3 x 10 1 N/m 2 

ami L = 47 in. The initial conditions used the assump- 
tion of isentropic recovery to obtain the cavity temper- 
ature while maintaining constant pressure across the 



Fig. 7: U.S. Army Airborne Optical Adjunct 

aperture. A characteristic inflow condition was used for 
the lip-blowing boundary, the flow rate computed using 
flight data and assuming isentropic compression of the 
ram air utilized in the aircraft. The 100% lip-blowing 
rate case corresponded to a iii = 0A2(pu) x . For the 
discussion below, computed high and low lip-blowing 
rates refer to 100% and 1% of this mass flow rate. The 
coarsened near-ficld grids arc shown in Fig. 8. 



Fig. 8: 2-d treated cavity: near field grids 

Tiie mechanism by which an aft ramp reduces the 
cavity feedback was explained by Heller and Bliss. 20 
Physically, this result can be explained from a force bal- 
ance normal to a streamline approaching the stagnation 
point. About the stagnation point, the velocity gradi- 
ent across the impinging shear layer creates a pressure 
gradient. However, there is a counteracting pressure 
gradient, pq 2 /r , due to the differing radii of curvature 
above and below the dividing streamline. 

For a rectangular cavity, the extreme of normal im- 
pingement of the shear layer onto the aft bulkhead 
causes further deflection into the cavity. Mass inges- 
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t ion into t lie cavity causes increased pressure, deflect- 
ing the shear layer outwards. With the shear region 
now outwardly deflected, mass expulsion from the cav- 
ity reduces the cavity pressure, inducing another cycle. 
Therefore, between the extremes of a normal or tan- 
gential impingement, of the shear layer, a balance of 
forces may be found. Use of a ramp instead of a con- 
vex surface at the reattachment region prevents shear 
layer perturbations from inducing instabilities of the 
type seen in rectangular cutouts. The length of the 
ramp must be large enough to accommodate the mag- 
nitude of the transverse shear layer excursions expected 
during operation. 

Computed and flight mean Mach number profiles are 
compared in Fig. 9 for two lip-blowing rates. The cpian- 
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Fig. 9: 2-d treated cavity: instantaneous Mach num- 
ber contours and mean profiles at (a) low and (b) high 
lip-blowing rate 

tity 0 indicates the angle from the cupola crest at which 
the data was measured. Figure 9 also shows Mach num- 
ber contours for the two lip-blowing rates above each 
set of profiles. The Mach number contours arc instan- 
taneous while the profiles were averaged over 2000 time 
steps. The difference between experiment and compu- 
tational results on the lower edge of the shear layer 
(r < -2") may be due to the computational simplifi- 


cation of an empty cavity. The difference at the upper 
aft portion of the shear layer ( r > 2”) appears to be 
due to blockage in the computational model. Overall, 
the maximum vorticity as a function of x-station is in 
agreement for both cases. 

Comparisons of low lip-blowing rate power spectra 
at the aft ramp are shown in Fig. 10. The computed 



Frequency, f, Hz 

Fig. 10: 2-d treated cavity, low lip-blowing rate: power 
spectra 

result generally falls 20 dB below the flight data, and 
a peak in the low lip-blowing rate spectra at 340 Hz is 
clearly computed, but is not seen in the flight data. 
Figure 10 also depicts data, the ordinate scaled by 
an( j the abcissa by * obtained from 

an AOA wind tunnel test. 21 The data can be seen to 
agree more closely with the computed results than with 
flight data, and a small peak exists at the computed 
frequency peak. The computed high lip-blowing spec- 
tra was broadband, but generally underpredicted flight 
data. 19 Computed power spectra were found using 4090 
points and a square window. 

Experimental data 13, 22> 23, 24 indicates that the fre- 
quency of large structures in shear layers occurs at 
Strouhal numbers of St — — = 0.024 ±0.003, where 

li| tUj 

9 is the local shear layer momentum thickness and / dc- 
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notes frequency. The computed Strouhal number can 
be estimated using Gortler's solution, given in Fig. 2a, 
to obtain 0 = 0.030]— ~r for rr 0 — 11.0, which is com- 
parable to the empirically determined correlation 2 ' of 
0 = 0.034j-”.r. Specifying r = 0.2 in this relation- 
ship along with a compressibility correction, 20 the com- 
puted peak in the AOA solution at 340 Hz corresponds 
to a Strouhal number of 0.032. Note that by using 
Rossiter's formula, given in Fig. 5, the frequency ob- 
tained for m = 5 is 3G0 Hz. From Fig. 11 it can be seen 
that four vortical cycles exist (m„ = 4), implying that 
— m — m v — 1. 
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Fig. 11: 2-d treated cavity, low lip-blowing rate: instan- 
taneous contours of 

Pc o 

The cause of the spike in the power spectrum was 
found from inspection of density fluctuations, which 
revealed large convecting shear layer structures. Fig- 
ure 11 shows a contour plot of depicting the growth 
and propagation of the shear layer structures. The 
large structures, associated with a 0.03 vertical ve- 
locity component, are the primary contributors of the 
computed density fluctuations of the shear layer. The 
speed of the structures is 0.56«oo» 19 below the value of 
O.GGu^o inferred by Rossiter for rectangular cutouts, yet 
above the O.oIhqo determined analytically by Roscoe 
and Hankev. 28 

Based upon these observations, it appears that large 
scale shear layer structures are being resolved. How- 
ever, the lack of empirical support from the flight data 
pressure power spectra is at odds with this conclusion. 
Three-dimensional effects are a possible explanation for 
the discrepancy. Rockwell 27 noted that for sufficiently 
large Reynolds numbers three-dimensionality reduces 
coherence in the shear layer. This implies that assump- 
tion of two-dimensionality for small flow oscillations 
may be suspect. In this two-dimensional computation 
the mass removed from the cavity by the shear layer 
entrainment process can only be replenished at the im- 
pingement region. This is in contrast to the mass ad- 
dition mechanism present in three-dimensions, which 
also includes mass replenishment via spanwise struc- 
tures such as streamwise vortices. The evolution of 


streamwisc-oriented vortirity interacting with the pri- 
mary vortices would act to spread peaks in the reat- 
tachment ramp pressure spectra. 

THREE-DIMENSIONAL RESONATING CAVITY 

An experimental investigation of SOFIA cavity qui- 
eting treatments was performed in the NASA Ames 
14' x 14' wind tunnel in 1990. 7 The geometry shown 
in Fig. 12 was the initial cavity configuration tested in 
the wind tunnel. This simulation was implemented in 
order to demonstrate the capture of self-excited cav- 
ity resonance in three-dimensions. The simulated flow 
conditions were the same as used in the wind tunnel, 
with the flow field initialized from a steady dean case. 
A stability-limited time step of At = 3.53 // s was used, 
corresponding to a streamwise CFL % 1 in the shear 
layer, and a CF L\„ iax « 500. 



Fig. 12: Configuration 25: wind tunnel and numerical 
models 

Sample pressure histories on the cavity walls and a 
comparison of the PSD resulting from the wind tunnel 
and numerical efforts are shown in Fig. 13. The PSD 
was obtained using 2048 points and a Hanning win- 
dow to match the treatment of the experimental data. 
The predicted frequencies of the dominant tones appear 
reasonable, and the computed dominant tone is within 
3 dB of experiment. The magnitudes of t ho computed 
higher modes are much lower t han observed experimen- 
tally. 







Fig. 14: Comparison of sound pressure levels 

if the same level of quieting could he predicted numer- 
ically as was observed experiment ally. Figure 15 shows 
the topology of the blended aft ramp grids used in this 
simulation of configuration 100. 


Fig. 13: Configuration 25: pressure histories and power 
spectra 

Estimation of the grid resolution required to main- 
tain a propagating wave of a specific magnitude can he 
deduced from the rectangular two-dimensional cavity 
and configuration 25 results. First, wavelength can he 
estimated by assuming the wave to he harmonic at a 
given frequency and travelling at the stagnation speed 
of sound. Next, it is noted that frequencies around 2 
kHz were resolved well in the two-dimensional case, in 
which the grid resolution was such that about 40 points 
supported the wave. From the configuration 25 results, 
it is seen that only the 700 Hz peak is well resolved, 
which again gives approximately 40 points across the 
wave for this coarser grid. Although numerical damp- 
ing of the higher frequencies can he expected, most of 
the energy is contained in the lowest frequency mode, 
as can he seen in the sound pressure level, or SPL , 
comparison of Fig. 14. Figure 14 also shows a quieted 
configuration which is discussed below. 


THREE-DIMENSIONAL QUIETED CAVITY 

During the wind tunnel experiments, use of a blended 
aft ramp resulted in the lowest sound production levels. 
This simulation was implemented in order to determine 


Telescope Primary^ 

cl 4^ 



Fig. 15: Configuration 100: cavity region topology 


For the SOFIA experiment, this type of geometry 
treatment was found to he quieter than the untreated 
configuration 25 case by over 30 dB. Figure 14 sum- 
marizes that t he proper trends were computed. The 
flow conditions were again initialized from a clean case, 
and integrated using a stability-limited time step size 
of Af = 7.0G//5. The frequency domain analysis was 
obtained using 409G points, a Hanning window, and no 
zero-padding. 
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A qualitative comparison of results ran he made from 
inspection of the computed instantaneous streamlines 
depicted in Fig. 14 and the oil flow from a similar 
molded geometry, shown in Fig. 16. A strong cross flow 
component is evident at the aft ramp for this aperture 
elevation angle, and is seen computationally as well. 



Fig. 16: Experimentally' observed surface flow pattern 

Quantitative comparisons were made for this pas- 
sively quieted cavity geometry in terms of shear layer 
profiles and pressure spectra in the cavity. Since errors 
in shear layer mass entrainment rate would adversely 
affect the cavity velocity field and hence the mean 
telescope loads, an important validation parameter is 
the shear layer spread rate. Figure 17 depicts mean 
experimental and computational shear layer Mach mini- 



Mach number 


Fig. 17: Configuration 100: instantaneous Mach number 
contours and mean profiles 

her profiles. The vertical scale of the profiles is twice 
that shown in the contour plot for clarity. Note that the 
experimental profiles were obtained using a rake, which 


is sensitive only to the x-rompnnent of velocity. Fig- 
ure' 17 indicates reasonable agreement in profile shape, 
however the computed growth rate is lower than seen 
in the data. This discrepancy may be due to geometry 
modifications in the experiment, or specification of an 
ovcrly-large value of (7q. 

Assessment of the oscillating telescope assembly loads 
requires the accurate resolution of the unsteady pres- 
sure field in the cavity. Comparison of the computed 
and observed spectra at specific locations provides a 
measure of confidence for the computed telescope loads. 
Toward the estimation of loads. Fig. 18 shows the pres- 
sure history and resultant PSD on the cavity walls. Al- 
though the peak levels are in agreement, the computed 
spectra can again be seen to drop more rapidly with fre- 
quency than the experimental results due to numerical 
dissipation. The spectral peak at 1800 Hz corresponds 
to a Strouhal number of 0.027, using a velocity ratio of 
0 . 1 . 



Fig. 18: Configuration 100: pressure history and power 
spectra on cavity wall 

THREE-DIMENSIONAL QUIETED AFT CAVITY 

Recently a new configuration which places the cavity 
aft of the wing, as shown in Figs. 19 and 20, has been 
proposed as a cost-effective configuration. However the 
cavity response of this aft telescope installation is of 
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Fig. 19: Aft telescope configuration: modelled geometry 

concern. Limited wind tunnel access has prevented 
timely resolution of this safety issue through exper- 
imental testing, hence this computation will provide 
an early measure of performance. Flight conditions at 
41,000 feet and M 0 0 = 0.85 were used, resulting in a 
Be i — 2.3 x 10'. Powered engines were simulated by 
specification of characteristic conditions for the inlets 
and both fan and core exhaust nozzles. The tail geom- 
etrv was not available for this simulation. 
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Fig. 20: Aft telescope configuration: cavity grids 

As with the forward cavity computations, a clean air- 
craft flow field was used as an initial condition for the 
cavity computation. Figure 21 compares the computed 
pressure coefficient profiles along the crest, side, and 
bottom of the model with flight data, 29 where the cav- 



ity location is shown for reference. Computed pressures 
are generally in agreement with the data, but the com- 
puted gradient, aft of the cavity location is more adverse 
than seen in the flight data. This is probably due to dis- 
crepancies in the empennage geometry and the lack of 
tail stabilizers. 
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Fig. 21: Clean configuration: pressure coefficient compar- 


Figure 22 shows the computed PSD on the aft rami) 
as compared to the scaled forward cavity wind tun- 
nel data. The spectra was computed using 1021 points 
stored at 5 At ^ 0.G8 ms intervals. The computed and 
scaled experimental spectra can be seen to be similar in 
shape and magnitude to about 100 Hz, where numerical 
dissipation becomes significant. The cavity grid cells. 






130 T 



Fig. 22: Aft telescope configuration: power spectra 


approximately 3 in, on a side, had been sized so that 
only low frequencies were well resolved. Use of smaller 
cells would capture short wavelength, high frequency 
waves which are unable to excite the massive telescope 
assembly. A Strouhal number of 0.028 is computed for 
the spectral peak seen at 110 Hz. 

Simulation of the distortion of a planar wave train 
of light propagating through the cavity aperture 19 re- 
veals the presence of convecting structures, shown in 
Fig. 23. The wavefront distortion, or optical path dif- 
ference (OPD), is proportional to the integrated den- 
sity normal to the aperture. Using f v — ^ and Fig. 23, 
it can be verified that the St = 0.028 peak is caused by 
the resolved shear layer structures. The structures con- 
vect with the local downwash velocity field caused by 
the wing. A large spanwise variation of the wavefront 
distortion can also be seen, with maximum < OPD ' > 
values at the crossflow edges of the aperture. Since the 
shear layer grid is uniformly distributed in this region, 
this spanwise variation may be indicative of increased 
vortex coherence at the boundaries of the aperture. 

CONCLUSIONS 

The development of an aerodynamic configuration 
analysis method has been summarized. The topology 
used within the overset framework has allowed simple 
and fast mesh generation while maintaining high reso- 
lution in critical flow regions. 

Evaluation of the overset mesh method was per- 
formed through comparison against experiment for sev- 
eral related flows. Sound pressure levels for resonant 
and quieted cavities compared well with experiment. 
However, use of a planar flow assumption for a quieted 
cavity appears to have caused a computed spectral peak 
where flight data indicated broadband behavior. Com- 
puted shear layer structures were found to cause the 



Fig. 23: Aft telescope configuration: wavefront distortion 

spectral peak at a Strouhal number of approximately 
0.03. In three-dimensions, a similar peak was seen in 
both experiment and simulation. Additionally, large 
spanwise density variation of the aperture structures 
was found. 
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